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Investigation  of  dynamic  fracture  of  long  brittle  beams  under  different 
loading  conditions  was  effected  under  ARO  Grant  DAAG29-78-G-0086.  The 
principal  investigator  was  Professor  G.  Herrmann  of  Stanford  University. 

Other  scientific  personnel  involved  in  the  research  were  Professor  M.  Shmuely, 
visiting  professor  from  the  Israel  Institute  of  Technology  (Technion)  and  Hr. 

C.  Levy,  graduate  student  at  Stanford. 

The  work  has  proceeded  in  two  distinct  areas  -  the  numerical  solution  of 
crack  propagation  problems  using  a  finite  difference  code  and  the  analytical 
solution  of  some  beam  fracture  problems  using  several  one-dimensional  theory 
formulations. 

The  motion  of  the  crack  tip  for  the  double  cantilever  beam  (DCB)  speci¬ 
men  and  for  the  single  edge-notch  specimen  (SEN)  was  obtained  for  several  load¬ 
ing  conditions  by  means  of  the  SMF2D  Code  proposed  by  Shmuely  and  Perl  In  a 
recent  paper.  The  numerical  results  for  the  DCB  model  were  compared  with  the 
experimental  results  of  Kalthoff  et  al.  and  found  to  be  within  good  agreement 
numerically  as  well  as  qualitatively.  Because  no  experimental  results  existed 
for  the  acceleration  stage  of  the  propagating  crack,  several  test  functions 
for  the  acceleration  stage  were  formulated.  These  functions  showed  no  appre¬ 
ciable  effect  on  the  acceleration  of  the  crack,  but  did  show  definite  effects 
on  the  deceleration  stage  of  the  propagating  crack. 

The  code  was  also  applied  to  the  SEN  model  for  which  experimental  results 
do  not  exist.  It  was  found  that  in  all  cases  the  dynamic  stress  Intensity 
factor  characteristics  depended  strongly  on  initial  crack  length  and  final 
crack  velocity.  All  cases  showed  the  characteristic  decrease  in  dynamic  stress 
Intensity  factor  as  the  crack  started  to  propagate,  thereby  confirming  the 
results  of  Freund  for  a  half-plane  crack  propagating  in  an  unbounded  medium 
(where  reflection  of  stress  waves  from  the  boundary  to  the  crack  tip  do  not 
play  an  Important  role). 
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Based  on  the  results  obtained,  the  SMF2D  Code  was  found  to  be  a  very  good 
numerical  tool  for  predicting  crack  behavior  within  the  limitations  of  the 
code.  As  a  result  of  this,  a  paper  was  submitted  for  publication  and  the 
work  on  the  numerical  solution  of  crack  propagation  problems  (M.  Shmuely  and 
C.  Levy,  "The  Dependence  of  the  Dynamic  Material  Toughness  on  the  Velocity 
and  on  some  Loading  Parameters")  is  now  in  print  in  the  international  Journal 
of  Fracture. 

The  second  area  of  endeavor  is  an  attempt  to  use  a  modified  quasi -static 
approach  to  the  fracture  of  infinitely  long  beams  under  the  assumption  that 
e  crack  propagation  velocity  is  much  smaller  than  any  of  the  characteristic 
&  speeds  of  the  beam.  The  problems  that  have  been  solved  analytically  have 
been  motivated  by  earlier  work. 

The  analytical  Investigation  of  dynamic  fracture  of  long  brittle  beams 
under  two  different  loading  systems  was  undertaken  assuming  that  the  crack 
propagation  velocity  was  smaller  than  any  of  the  characteristic  wave  speeds 
of  the  beam.  This  assumption  allowed  us  to  formulate  the  problems  via  a 
modified  quasistatic  approach.  To  make  the  problems  more  tractable,  the  one¬ 
dimensional  beam  models  of  Euler-Bernoul 1 i ,  Rayleigh  and  Timoshenko  (wtthout 
rotary  Inertia)  were  incorporated  in  the  problem  formulation.  The  loading 
(statically  applied)  was  either  remote  axial  tension  or  remote  pure  bending 
and  the  results  obtained  were  the  crack  tip  velocity-time  history;  the  crack 
tip  displacement-time  history;  the  bending  moment  and  the  extensional  load  on 
the  fracturing  cross-section.  To  preclude  reflection  effects,  the  beam  was 
considered  to  be  Infinite  in  length.  The  results  obtained  showed  that  frac¬ 
ture  is  not  complete  for  either  applied  loading,  but  requires  the  accounting 
of  the  reflection  process  for  fracture  to  be  total.  The  following  results  do 
not  account  for  the  reflection  process,  but  this  feature  Is  now  under  Investigation. 
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The  simplest  beam  model  studied  was  that  of  Euler-Bernoull l .  The  analyt¬ 
ical  formulation  for  the  beam  under  remote  applied  bending  moments  may  be 
found  in  the  literature.  However,  an  error  ex  sts  in  the  formulation  of  the 
problem  and  when  corrected,  brings  th;  results  closer  to  earlier  experimental 
observation.  To  investigate  the  rotary  inertia  effects  on  the  bending  moment 
and  extenslonai  loading  at  the  fracturing  cross-section,  the  Rayleigh  beam 
model  was  employed.  With  respect  to  the  results  of  the  Euler-Bernoull I  model, 
for  a  given  initial  crack  length,  the  rotary  inertia  tended  to  retard  crack 
growth,  to  decrease  the  bending  moment  and  to  increase  the  extenslonai  load 
on  the  fracturing  cross-section;  and  to  reduce  the  crack  tip  velocity  over 
the  majority  of  the  beam  height. 

The  final  model  Investigated  was  the  Timoshenko  beam  model.  The  Timo¬ 
shenko  model  was  employed  in  order  to  show  the  shear  effects  on  the  data  to  be 
obtained.  With  respect  to  the  results  of  the  Euler-Bernoul 1 1  model  for  a 
given  Initial  crack  length,  the  shear  tended  to  enhance  crack  growth;  to  in¬ 
crease  the  bending  moment  and  to  decrease  the  extenslonai  load  on  the  fractur¬ 
ing  cross-section;  and  to  Increase  the  crack  tip  velocity  over  most  of  the 
beam  height.  Thus  by  "averaging"  the  data  of  these  last  two  models  we  found 
that  the  overall  results  were  either  close  to  the  Euler-Bernoul II  model  or 
favored  very  slightly,  the  Rayleigh  beam  model  data.  This  "averaging"  could 
be  considered  a  first  order  approximation  to  the  Timoshenko  beam  model  and 
the  results  tended  to  show  that  the  Euler-Bernoul 1 i  theory  provided  us  with 
an  adequate  model  for  the  dynamic  fracture  of  infinitely  long  notched  brittle 
beams  under  varied  loading  conditions,  in  all  It  is  found  that  all  the  models 
adequately  describe  the  fracture  process  unti 1  wave  reflection  effect  came  into 
play.  This  second  area  has  resulted  in  six  manuscripts  that  have  been  sub¬ 
mitted  for  publication: 
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On  the  Effect  of  Axial  Force  on  Dynamic  Fracture  of  a  Beam  or  Plate  In 
Pure  Bending 

Dynamic  Fracture  of  a  Beam  or  Plate  under  Tensile  Loading. 

Symmetric  Dynamic  Fracture  of  a  Beam  or  Plate  under  Tensile  Loading 
Effect  of  Bending  Moment  on  the  Dynamic  Fracture  of  a  Beam  or  Plate  under 
Tens! le  Loading 

Effect  of  Shear  and  Rotary  Inertia  on  Dynamic  Fracture  of  a  Beam  or 
Plate  under  Tensile  Loading 

Effect  of  Shear  and  Rotary  Inertia  on  Dynamic  Fracture  of  a  Beam  or  Plate 
In  Pure  Bending 
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